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ABSTRACT. We address several issues regarding the interpretation of galactic 7 -ray sources. 
We consider powerful pulsars in binaries producing X-ray and gamma-ray unpulsed emission 
from the shock interaction of relativistic pulsar winds with circumbinary material. Nebular 
mass outflows from companion stars of binary pulsars can provide the right calorimeters 
to transform a fraction of the electromagnetic and kinetic energy of pulsar winds into high 
energy radiation. We discuss the physics of interaction of relativistic pulsar winds with 
gaseous material and show that the conditions in pulsar binary systems might be ideal to 
constrain shock acceleration mechanisms and pulsar wind composition and structure. We 
briefly discuss the example of the 47 ms pulsar PSR 1259-63 orbiting around a massive 
Be star companion and monitored by X-ray and gamma-ray instruments during its recent 
periastron passage. In addition to young pulsars in massive binaries, also a class of recycled 
millisecond pulsars in low-mass binaries can be interesting high energy emitters. 

1. Introduction 

Recent data of the EGRET instrument on board of the Compton Gamma-Ray 
Observatory (CGRO) revealed at least ~40 point-like 7 -ray sources of low galactic 
latitude bn ~ 10° (Fichtel et ah, 1994). EGRET conhrms several of the previously 
discovered COS-B sources (Swanenburg et ah, 1981; Bignami & Hermsen, 1983) 
and adds a signihcant number of new ones. Gamma-ray emission (from ~MeV to 
tens of GeV) is a manifestation of particle acceleration and non-thermal behavior of 
astrophysical systems. Candidate systems for high energy emission include: 

( 1 ) isolated pulsars, 

( 2 ) supernova remnants, 

(3) pulsars moving in the interstellar medium; 

(4) compact stars (e.g., pulsars) with relativistic winds in binaries, 

(5) massive binaries with colliding non-relativistic mass outflows, 

( 6 ) accreting compact stars. 
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Revealing the nature of unidentified gamma-ray sources is difficult, and requires 
monitoring and multiwavelength observations. For example, the 7 -ray source 2 CG 
195+04 (Geminga) (e.g., Bignami & Hermsen, 1983), was recently identified as a 
237 ms 7 -ray pulsar (Bertsch et al., 1992) only after a crucial ROSAT observation 
(Halpern & Holt, 1992). Accreting sources are not strong 7 -ray emitters (e.g., the 
black hole candidate Gyg X -3 is marginally detected by GOMPTEL up to ~10 MeV, 
e.g., Hermsen, 1994), a fact related to the difficulty of sustaining non-thermal acceler¬ 
ation and radiative processes in accretion disks. The environment surrounding com¬ 
pact stars must therefore be relatively ’clean’ to allow the formation and propagation 
of 7 -rays. Without invoking the existence of black hole related phenomena, pulsar 
magnetospheres (of 7 -ray pulsars) and nebular shock environments are systems con¬ 
sidered in a ’conservative’ approach to the interpretation of galactic 7 -ray sources. We 
mainly discuss here high energy shock emission in binary systems containing rapidly 
rotating pulsars interacting with circumbinary material. The high energy emission 
might be modulated with the underlying orbital period. These systems (both low- 
mass and high-mass binaries) are interesting candidates for the low-galactic latitude 
sources of Fichtel et a/.(1994). 

We address a few topics related to the interpretation of galactic 7 -ray sources in 
term of a pulsar-driven model of relativistic nebular emission. We do not discuss here 
isolated 7 -ray pulsars (which are an important class of steady 7 -ray sources, e.g., 
Thompson, 1994), pulsars moving in the interstellar medium (e.g., PSR 1957+20, 
Fruchter et al.^ 1988; Arons & Tavani, 1993, hereafter AT93), OB associations (e.g., 
Montmerle, 1979), and massive binaries with colliding non-relativistic mass outflows 
(e.g., the Wolf-Rayet system WR I4O, Williams et al.^ 1990; Fichler & Usov, 1993). 
The physics of interaction of relativistic winds with nebular material is of general 
interest, and it may also be relevant for transient high energy emission of 7 -ray 
bursts (Meszaros & Rees, 1992) and soft gamma-ray repeaters (Tavani, 1994d). 

2. High Energy Emission from Plerions: the Crab Nebula 

Pulsars lose their spindown energy mainly by the emission of electromagnetic 
waves and kinetic energy of relativistic particles. In the absence of a relatively dense 
environment, these relativistic outflows are practically unobservable, as demonstrated 
by the lack of detectable high energy emission from the majority of pulsars. How¬ 
ever, in a few cases, either a supernova remnant, a dense interstellar medium, or a 
dense circumbinary material, can provide the right nebular environment to convert a 
fraction of pulsar spindown energy into observable high energy emission. The nebula 
acts as a calorimeter with conversion efficiency £„ into high energy emission in the 
range 0.01 ^ Ca ^ 0 . 1 . 

Until recently, only nebular environments of pulsars embedded in supernova rem¬ 
nants were available to study. The Crab Nebula is the prototypical example of an 
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X-ray and 7 -ray radiating pulsar-driven nebula (plerion). The Crab Nebula shows un¬ 
pulsed emission from radio to 7 -ray wavelengths^. The 33 ms Crab pulsar is embedded 
in a quasi-spherical ’pulsar cavity’ with its relativistic particle and electromagnetic 
wind conhned by a MHD shock at the distance^ from the pulsar ~ 3 -10^^ cm (Rees 
& Gunn, 1974; Kennel & Coroniti, 1984, hereafter KC84). Observations show that 
the Crab nebula component dominates (by a factor ~5) the Crab pulsar component 
of the spectrum up to a few MeV with photon index cti ~ 2.2 (e.g., Bartlett et 
al., 1994). The nebular becomes softer (and subdominant with respect to the pulsar 
emission) with a 2 ~ 2.9 between a few MeV and a few hundreds of MeV (Nolan et 
a/.1993). At ~500 MeV a change in spectral index was observed by EGRET up to 
~30 GeV with ~ 2.7 (Nolan et al.^ 1993, De Jager et al.^ 1994). The photon index 
of the Crab nebula at the highest EGRET energies is consistent with the observed 
level of unpulsed TeV emission (Vacanti et al.^ 1991). 

High energy emission of the Crab Nebula is a manifestation of a relativistic pulsar 
wind originating from the pulsar (Rees & Gunn, 1974). As the pulsar wind (a mixture 
of electromagnetic helds and particles) interacts with a quasi-static MHD shock at r^, 
a redistribution in momentum space of the downstream particle distribution occurs 
and non-thermal synchrotron and inverse Compton (IC) radiation can be emitted'^. 
The volume of the emitting region is larger than the inner pulsar cavity, and both 
particles and helds diffuse out from the central region into the main body of the 
nebula. The Crab Nebula is the most thoroughly studied plerion and any model 
of the Crab nebular shock emission is a natural candidate for the study of pulsar 
shock emission of other systems. In this section, we use the Crab nebular model to 
introduce the main relevant properties of pulsar winds and nebular emission that we 
will apply below to binary systems. 

2.1 RELATIVISTIC PULSAR WINDS 

Pulsar magnetospheres are powerful particle accelerators and a rapidly rotating 
pulsars as the Crab is capable to emit a powerful wind from the open held line regions. 
The accelerating voltage drop above a neutron star polar cap region can be very high. 
For the Crab pulsar is ~ VE p/G ~ Volt, with 0 = 2 7 r/P the spin 

frequency of a pulsar of period P, and p = /iso 10^° (in c.g.s. units) the neutron star 
magnetic moment. Acceleration and e^-pair creation inside the light cylinder radius 
[Ric = c/0 ~ 1.57-10® cm for the Crab pulsar) leads to the formation of a relativistic 
particle wind. At sufhciently large distance from the pulsar, the magnetic held Pi is 
mostly toroidal and the MHD approximation can be applied to the pulsar wind. The 

^For a summary, see Kennel & Coroniti (1984) and references therein. 

^The MHD inner shock radius of the Crab Nebula has an apparent size of 10" at the distance of 
2 kpc. For all practical purposes, pulsar cavities of known systems appear as point sources to X-ray 
and 7 -ray instruments. 

^See however Cheung & Cheng (1993) for an alternative model. 


4 



wind has a flow rate in ions Ni and e^-pairs N±^ a flow Lorentz factor 71 , and a ratio 
of the electromagnetic energy flux (Poynting flux) to kinetic energy flux upstream of 
the shock, a = + m±(A^_l_ + A^_)] 7 ic^. A MHD model of a pulsar wind 

can therefore be characterized by a set of five parameters 

(cr,7i,5i,p,/p±,Z/A) (1) 

where we defined the ions’ density py the e^-pairs’ density and the charge to 
mass ratio for ions ZjA. 

The details of acceleration, particle composition and electromagnetic structure 
of relativistic winds near Ric are poorly constrained at present. Qualitatively, the 
quantity a is expected to be relatively large near the pulsar, i.e., a ~ 10 “^, with a rate 
of pair creation N± ~ 10^® s“^ and Lorentz factor ~ 200 (Coroniti, 1990). As the 
wind expands out from the pulsar, it evolves from a high-cr to a low-cr configuration 
and the 7 -factor increases. No direct information is available to probe the inner region 
of a pulsar wind near Ric- However, a region very far from Ric (with VsjRic ~ 10® —10® 
in the case of the Crab Nebula) can be (indirectly) probed by the properties of a MHD 
shock. An important result is obtained by combining the observational results of the 
Crab nebula with a detailed hydrodynamic and radiative model of the shock. With 
the assumption of a downstream power-law e^-pair momentum distribution, the MHD 
model of the pulsar wind in the Crab Nebula gives a ~ 0.005 and 71 10® (KC84). 

Notice that these values are characteristic of the Crab pulsar wind only for a very 
large distance from Ric- The determination of the ions’ content of the Crab pulsar 
wind is more controversial and dependent on the shock acceleration mechanism used 
to model the Crab perpendicular shock‘d. A particle-in-cell numerical simulation 
of a perpendicular shock for a pure e^-pair pulsar wind fails to produce a non- 
thermal downstream particle distribution (Hoshino et al.^ 1992, hereafter HAGL92; 
Gallant & Arons, 1994). A non-thermal downstream momentum distribution of e^- 
pairs develops only with the inclusion of a heavy ion component of the pulsar wind 
dominating the particle kinetic energy flux (HAGL92). Despite the one-dimensional 
character of these simulations and the approximations of the calculations of HAGL92, 
the possible presence of ions in the pulsar wind should be considered as an interesting 
possibility that can be tested by observations. 

2.2 NEBULAR SHOCK EMISSION 

The high energy emission produced by the pulsar wind termination shock of the 
Crab Nebula depends on shock acceleration mechanism. The fundamental fact about 
Crab-like pulsar-driven shocks is that the geometry of the shock for a spherically sym¬ 
metric outflow is typical of a ‘perpendicular shock’ rather than of a ’parallel shock’. 

■^A perpendicular shock is characterized by a shock surface whose normal vector is perpendicular 
to the upstream magnetic held, as expected in the idealized model of the Crab pulsar wind with a 
toroidal magnetic held conhned by a spherically symmetric inner cavity. 
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The particle acceleration mechanism is therefore different from the (hrst-order) dif¬ 
fusive Fermi acceleration mechanism that is applicable to ‘parallel shocks’ (e.g., Bell, 
1978; Blandford & Ostriker, 1978; Jones et al., 1994). Shock drift acceleration models 
of perpendicular shocks might explain pure e^-pairs acceleration in the Crab Nebula, 
but only if cross held diffusion is present (e.g., Jones et al.^ 1994). A developed 
model of perpendicular shocks currently available is the model of HAGL92, that we 
discuss here. In order to explain the properties of the Crab Nebula, it is necessary to 
identify a shock acceleration mechanism able to transform a relativistic Maxwellian 
particle distribution upstream of the shock into a downstream power-law distribution 
of e^-pairs. The downstream distribution of e^-pairs per energy interval becomes, 
N{'^)dq ~ 7 “®, with s an exponent. Simulations indicate that a resonant mech¬ 
anism of e^-pair acceleration mediated by ions’ gyromotion satishes the necessary 
requirements for the Crab Nebula (HAGL92). 

Once a downstream non-thermal particle distribution is established, radiation 
occurs by synchrotron emission and inverse Compton scattering. The magnetic held 
of the Crab Nebula upstream of the inner shock is 



Er 


1/2 

~ 3 • 10“® 


a 




-|l/2 


0.005 (1 + 17) 


G 


( 2 ) 


with fp the fraction of the total solid angle for the pulsar wind expansion [fp = 1 for 
isotropic expansion). In Eq. (2) we used a typical combination of density and gas ve¬ 
locity at the Crab Nebula shock, n 2 = n(rs)/(10^ cm“‘^) and v-j = n(rs)/(10^ cms“^). 
Note that Bi[rs) depends ultimately on the value of the gas pressure at and, to 
hrst order, is independent of the pulsar spindown power. In the case of the Crab 
Nebula, the nebular pulsar pressure is calculated by taking into account the energy 
accumulation during the pulsar lifetime (Rees & Gunn, 1974). The Rankine-Hugoniot 
equations of a relativistic shock give for a low-u pulsar wind the value of the down¬ 
stream magnetic held i ?2 = 3 i?i (KC84). The downstream magnetic held B 2 enters 
in the calculation of the shock synchrotron emissivity. 

The observed efhciency of conversion of pulsar spindown energy into nebular high 
energy emission is near ~10% (KC84), and the numerical model of HAGL92 is in 
agreement with this relatively high efhciency. The high energy spectrum of the Crab 
Nebula is explained in terms of synchrotron radiating e^-pairs being accelerated at 
the inner shock and progressively diffusing out in the central nebula (KC84). Inverse 
Compton scattering against a plerionic diffuse low-energy photon background may 
play a role for the nebular emission above ~500 MeV (De Jager et al.^ 1994). 


2.2 A CA+IANAAN EXPERIMENT: COMPRESSING PULSAR-DRIVEN CAVITIES 


A pulsar cavity in a supernova remnant is by dehnition quasi-static without ap¬ 
preciable time variability. If pulsars in supernova remnants were the only systems 
available to study pulsar winds, there would be no hope to see how a MHD shock 
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changes its radiative properties as the pulsar cavity changes in size. However, one 
can think of a gedanken experiment with pulsar cavities. Let us imagine a rapidly 
rotating and energetic pulsar at the center of a cavity in a plerionic and calorimetric 
nebula. The nebular high energy emission depends on the pulsar wind parameters of 
Eq. (1). Given the relatively large number of (model-dependent) parameters, it might 
be difficult to obtain a unique model of emission. Now, let us imagine to compress 
progressively the pulsar cavity (e.g, by increasing the gas pressure at the MHD shock) 
from a large radius with VslRic ~ 10® as in case of the Crab Nebula to a much smaller 
radius. The local magnetic held strength varies, and the corresponding change of the 
high energy intensity and spectrum will follow a pattern proper to a particular set of 
pulsar wind parameters. For example, the magnetosonic e^-pair acceleration mecha¬ 
nism for a perpendicular shock mediated by heavy ions can be tested for comparable 
or less® than the ions’ Larmor radius, ~ 3 • 10®^ (A/Z) (7i/10®)(i?i/10“® G)“®. 
A transition to a power-law high energy emission to a quasi-thermal spectrum of 
emission is expected in the model of HAGL92 as the shock radius decreases and the 
magnetosonic mechanism becomes progressively more inefficient for R^g ~ r^. Anal¬ 
ogously, a transition from a very-low-cr to a moderately-low-cr pulsar wind might in 
principle produce observable effects, as decreases. Therefore, time dependent nebu¬ 
lar environments can provide ideal calorimeters to test theories of shock acceleration 
and pulsar wind composition. 

Obviously, it is not possible to perform this gedanken experiment with plerionic 
supernova remnants. However, pulsars in binary systems can be in the right environ¬ 
ment to produce time variable nebular emission. 

3. Pulsar-driven Nebular Emission of Binary Pulsars 

Binary pulsars are interesting systems which may provide environments where we 
can test the pulsar wind model and acceleration mechanism established for the Crab 
Nebula. There is no reason in principle to assume that the pulsar wind parameters 
of Eq. (1) are the same for all pulsars. In particular, binary pulsars are characterized 
by a ratio VslRic different from the Crab value. 

3.1 PULSAR-DRIVEN SHOCKS IN CIRCUMBINARY MATERIAL 

We consider here a binary containing a pulsar interacting with a mass outflow 
(gaseous non-relativistic wind from the companion star surface). The discussion 
presented here can be applied to high-mass and low-mass systems with proper mod- 
ihcations. We follow here the discussion of Tavani, Arons & Kaspi (1994; hereafter 
TAK94). 

The mass outflow has in general a density prohle p{R) = po{RlR*)~^ , where i?* is 

®We notice that for a radial dependence of electromagnetic pulsar waves Bi ~ the ratio 

Vg/RL,i varies only because of a change of a with radius. 
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the radius of the star, R the distance from its surface, and po = Af/[4 7 r/i?^ Uo], with 
M the mass outflow rate, Vq the outflow velocity at the surface of the companion star, 
and / the fraction of dvr steradians into which the mass loss occurs. Typically, the 
surface density po is of the order of 10“^^ —10“^^ g cm“^, for both massive Be-star-like 
companions and low-mass companions. For a constant velocity mass outflow (typical 
of low-mass star outflows), the ‘outflow exponent’ is n' = 2. More generally, for 
massive companions, the outflow exponent is in the range 2 < n' < 4 (Waters et al., 
1988). The mass outflow pressure at distance R from the surface of the companion 
star has a radial dependence of the form, P{R) = p{R) n(i?)^ {RjR^Y ^5 obtained 
with the use of the continuity equation in spherical symmetry. Note that the outflow 
velocity changes with distance as v{R) ~ {RjR^"' i-e., it increases for n' > 2 

as suggested by observations of Be mass outflows (Waters et al.^ 1988). Near the 
pulsar orbit, the outflow velocity can reach values of order of 100 kms“^. We dehne 
the mass outflow parameter®, T = Af-sCoy, where M = (10“® Af 0 yr“^) Af_ 8 , and 
Vq = (10® cms“^)no, 6 . It is useful to parametrize the properties of the companion 
star mass outflow by specifying the parameters (T, n'). 

The radial distance of the apex of the shock region from the pulsar is obtained 
by equating the pulsar radiation pressure with the ram pressure of the mass outflow 
(e.g., AT93). The radius is the solution of the equation (TAK94), 

rsIdjO) ^ f ^ 7.2f/^ 

[1 - (r,M 0 ))] 2 -W 2 [cMvqJ [d(0)J Ti /2 [d(0)J ’ 

where we have assumed i?* ^ d, with d the pulsar-companion star separation. In 
general, is a function of the true anomaly d, n', and T. Pulsar radiation pressure 
of a rapidly rotating pulsar can resist the compressing tendency of mass outflows in 
binaries and avoid accretion (e.g., Tavani, 1991). This conclusion follows from the 
fact that for high-mass binaries the gravitational accretion radius 

Rg = 2GM,lvl ^ (3 • 10'2 cm)M,,i.4/<7 (4) 

(with Afp^i .4 the pulsar mass in units of 1.4 Af 0 and = n*/(10^ cms“^), where 
n* = and Vorb is the orbital velocity) is smalP compared to the shock 

radius r^ for a large range of outflow parameters T and n'. However, the effects of 
radio pulse dispersion and absorption as the pulsar passes into the densest part of the 
outflow may suppress the pulsed radio signal. The pulsar can be totally or partially 
‘hidden’ or enshrouded during its orbital motion. 

The magnetic held Bi{rs) at the apex of a pulsar cavity in a typical binary is 
given by Eq. (2) with n 2 {rs) ~ 10® — 10®. 

®We use the notation Xp, for a physical quantity in units of 10^. 

^This conclusion is not true in general for low-mass binaries, and pulsar cavities of low-mass 
systems can be subject to interesting instabilities (Tavani, 1991). 
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4. Shock Radiative Regimes of Binary Pulsars 

The luminosity and spectrum of high energy shock radiation depend on the var¬ 
ious radiative timescales and their relation to the flow timescale. The flow time of 
particles streaming with velocity Vj ~ c/3 behind the nebular shock is Tf ~ — 

50 [l/(uio/p [('^/3)/h/] s. We adopt here the acceleration time of nonthermal 

magnetosonic acceleration of e^-pairs ~ 0 ”^ = rriiCji j{ZeB 2 ) ~ 10 ^ -^ 2 ] s, 

where Oj- is the ion’s relativistic cyclotron frequency at the shock front (HAGL92). 
We can assume the validity of the perpendicular shock calculations for the apex of the 
pulsar cavity. These calculations give a downstream power-law index of the e^-pair 
momentum distribution s ~ 2 , and a fraction £„ ~ 0.2 of the flow energy flux going 
into the hnal energy of the nonthermal pairs (HAGL92). However, for a variety of 
binary orbital parameters, because of the rapid cooling of the shocked relativistic 
particles, near periastron the downstream relativistic Lorentz factor never reaches 
the value 7^ = 71 m ±) which is obtainable only without appreciable radiative 

cooling. For a very eccentric orbit (as in the case of PSR 1259-63 (Johnston et al., 
1992) only far from periastron the downstream 7 -factor can reach 7 ^. 

Radiative losses of the e^-pairs due to synchrotron radiation in the local magnetic 
held at the shock and inverse Compton scattering (IC) in the background IR, optical 
and UV radiation held originating from the companion star surface may limit the 
maximum Lorentz factor attainable at the shock. This is particularly true for high- 
mass pulsar binaries, with luminous companion stars providing a large background 
of low-energy photons permeating the shock region. The synchrotron loss time is 
Ts( 7 ) = 3m^c®/(e‘^i?27) ~ 1-3 • 10^ [ 70 -^ 2 ]”^ s- The cooling time for IC scattering of 
high energy e^-pairs against photons in the Thomson regime (7 e ^ m± c^, with e the 
photon energy) is Tc = {Ub lUp)Ts^ where Ub is the energy density of the post-shock 
electromagnetic held, and Up is the radiation energy density at the shock radius. 
The electromagnetic energy density at is UB{rs) = i? 2 (rs)/ 87 r, and the radiation 
energy density is Up = R^^ctbT'^ j[R^c), where cr^ is the Stefan-Boltzmann constant, 
T is the companion star’s effective temperature, and Rs is the distance between 
the closest point on the surface of the companion star and the shock. The above 
formulae are valid only in the Thomson regime and in the case of a radiation bath 
at the photospheric temperature of the companion (e.g., T ~ 27,000 K ~ 2.3 eV for 
a typical Be star). This approximation is valid only for 71 ~ 10®. However, for a 
typical Crab-like Lorentz factor of the relativistic particles at the shock (71 ~ 10®), 
Compton scattering occurs in the relativistic regime and it is necessary to take into 
account the Klein-Nishina decrease of the cross section as a function of center-of-mass 
energy. The dimensionless parameter characterizing the scattering is (Blumenthal & 
Gould, 1970) Te = 47 e/(m± c^), with e = 2.7 Ub T the average photon energy of the 
radiation blackbody held, and kB the Boltzmann’s constant. For typical Crab-like 
parameters of the pulsar wind and Be-star surface temperatures, F^ ~ 50 71 ^ 6 - The 
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Thomson approximation can be used for T^ ~ 1, but it is clearly not applicable 
if 7 i,6 ~ 1- It can be shown that the ratio of the IC scattering timescale in the 
Klein-Nishina regime over the non-relativistic Thomson timescale Tc in the case 
of PSR B1259—63 with 71 ~ 10® is ~ 10. 

Different radiative regimes are produced at the shock for different combinations of 
timescales. A value of small compared with other relevant timescales is required to 
accelerate the e^-pairs and produce a non-thermal power-law distribution of at the 
shock with energies 7 > 71 . On the other hand, the cooling processes by synchrotron 
and IC emission limit the non-thermal shock acceleration. The total cooling rate 
per radiating e^-pair at the shock can be written as . If 

is long compared with tj, the radiative losses are negligible, the flow is adiabatic, 
and the resulting radiation nebula is diffuse (e.g., AT93). On the other hand, if the 
cooling timescale is short compared to tj, the shock acceleration of pairs produces a 
spectrum of which is cut off at high energy, and a compact nebula forms. Finally, 
if Ta ^ Tr, uo uouthermal acceleration is possible at the shock and the emerging 
spectrum is dominated by the cooling mechanism with the largest rate. 

If heavy ions constitute a substantial component of the pulsar wind, the max¬ 
imum Lorentz factor 7 ^ that can be achieved by nonthermal acceleration of ions’ 
magnetosonic waves at the shock is obtained from the condition Ts{rs)lTa{rs) ~ 
(hm 7 i)~^ which is substantially smaller than 7 ^ obtained in the absence of radiative 
losses. Thus, different radiation regimes occur depending on the ratios of timescales 
Ts/tc, Ta/vs, Ta/vc. We uotice that Tc can be much smaller than near periastron for a 
large variety of outflow parameters. IC scattering therefore can dominate the cooling 
processes at periastron unless the shock is relatively close to the pulsar. This hap¬ 
pens for T ~ 10® and typical high-mass binary pulsar orbital parameters (TAK94). 
A pulsar-driven shock radiative regime dominated by IC cooling due to a companion 
star optical emission has not been studied previously. 

4.1 INTENSITY AND SPECTRAL PROPERTIES OE SHOCK EMISSION 

4 . 1.1 Non-thermal ‘diffuse’ nebular emission. This is the case characterized by 

Ta Tf Nl Ts^Tc. This case is relevant for a fraction or whole pulsar orbits with 
orbital periods Porb ~ lyr. Except for cases with 1 r-^ r-^ 10 (with n' = 3, 

and cr 1 ), synchrotron cooling dominates the emission for negligible emission 
from the companion star or large orbital distances. Since the flow time in the 
shock front is much less than the synchrotron loss time in the downstream flow, 
the flow is adiabatic both in the ions and in the pairs. Therefore, the resulting 
synchrotron nebula at distances from the pulsar r ~ is diffuse. As a result, the 
shocked wind radiates a nebular synchrotron spectrum with emissivity jc oc 
(photon energy cm“® s“^). The nonthermal acceleration of pairs at the shock pro¬ 
duces a power-law distribution of accelerated particles for the case of magnetosonic 
absorption of the ion’s waves and the maximum Lorentz factor 7 ^. Therefore, 
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the energy spectrum is in the range Ci < e < e^, where Ci = 0 . 37 ^ha;c(i? 2 ) — 
(0.05 keV) {Z/Ay {ng (cr/0.005)^/^, with iOc the electron/positron cyclotron fre¬ 

quency, and Cm = j'yfiiOc{B 2 ) — (210 MeV) (us Cg)^/^ (cr/0.005)^/^ (TAK94). Since 
most of the emission comes from the highest energy particles, the nonthermal part of 
the pair spectrum dominates the emission. For iV( 7 ) oc 7 “®, the bolometric luminos¬ 
ity of the flowing plasma is sa (7 t{B2 jG tt) (1 — x^~^)/[l — 

where aj is the Thomson cross section, the number of pairs in the nonthermal 
part of the pair distribution function, Vg = {kr^y (with A; ~ 1 ) is the volume occu¬ 
pied by the radiating wind downstream from the shock, and x = = Zm±lmi 

(AT93). We estimate the total number of suprathermal pairs by noting that the 
shock converts a fraction £„ of the pulsar’s energy loss into the energy of the non¬ 
thermal part of the pair spectrum. Therefore, it can be shown that in the downstream 
medium flowing with velocity Vj, Lyj ~ 10^“^ (cr/0.005) (ea/0.2) (ns vH erg s“^ . 

4-1.2 Non-thermal ‘compact’ nebular emission. This is the case characterized by 
Ts, Tc ^ Tf which applies to tight binaries or to the pulsar passage near periastron for 
systems similar to PSR 1259-63. We distinguish two main radiative regimes (TAK94). 

Case A 

If Ta <C. Tf, Tg, Tc, uou-thermal acceleration of e^-pairs at the shock is possible and the 
downstream power-law distribution extends from 71 to some This case occurs only 
for a small Vg near the pulsar and for a large ram pressure from the companion star 
outflow, i.e., T ~ 100 for high-mass binaries. Synchrotron radiation may dominate 
the cooling near periastron with ^ Tg for T ~ 10^, especially if 71 ~ 10® when the 
relativistic IC scattering dominates. In this case, the energy spectrum is peaked at 
Cl ~ 0 . 37 ^ JiuJc{B 2 ) — 10 keV and it is a power law spectrum of emissivity ~ 
up to an energy ~ 40 MeV. 

Case B 

For values of T ~ 10®, IC scattering dominates the cooling in high-mass binaries® and 
the acceleration time becomes larger than Tg (and Tc,r)- In this case, the energy dis¬ 
tribution of downstream e^-pairs resembles a relativistic Maxwellian with no power- 
law tail due to shock acceleration. The optical depth r* for non-relativistic IC scat¬ 
tering in a high-mass system can be estimated as r*,, ~ ax Up Vg/e ~ 10 rgX 2 , where we 
used the estimate for the radiation energy density Up (40ergcm ®) [(R*/Rs)/0.2]®. 
If IC scattering is predominantly in the Thomson regime, the photon spectrum has a 
Maxwellian shape with most of the energy radiated near the peak energy Ep ~ 7 ®e. 
For example, in the case of PSR B1259-63, the center of mass energy in the electron 
frame is e 7 i/(m±c®) ~ 4. 571^6 and the Thomson regime applies when 71 ~ 10®. In 
this case, the peak energy is Ep ~ (45 GeV) 7 ^ 5 . The total luminosity of the IC spec¬ 
trum is Lc,nr = (4/3) (Tj cUp f N k where A; is a constant of order unity, N ~ Nf 

®IC cooling is typically negligible in low-mass binaries (e.g., AT93). 
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the total number of e^-pairs IC radiating per unit time, with N f En/^fi m± the 
e^-pair production rate, f the residence time f ~ rs/(c/3), and / the fraction of N e^- 
pairs near the shock. The total number of e^-pairs is ~ 10 ^®rs^i 2 Jil, and from an 
estimate of the solid angle of particles streaming toward the companion star, / ~ O.f. 
We obtain Lc',nr ~ 5 • 10^®(//0.1) Ws 71,5 ^ 3,12 erg s“^ . Alternately, if 71 ~ 10® as in 
the case of the Crab pulsar (KC84) the IC scattering is in the Klein-Nishina regime 
with maximum photon energy E^ax ~ 7 i m± ~ (500 GeV) 7 i, 6 - The IC lumi¬ 
nosity is in this case, Lc\r ~ (3/8) aj c Up {m± c^/e)^ / N [In (2 71 e/m± U) + 1 /2] ~ 
6 • 10®“^ [(i?*/i?s)/0.2]^ (//O.l) erg s“^ . Note that the e^-pairs of the pulsar wind can 
have their 71 lowered by IC scattering even before they reach the shock (TAK94). 

5. Applications to Pulsar Systems 

5.1 YOUNG PULSARS IN HIGH-MASS BINARIES: THE CASE OE PSR B1259-63 

As binary stellar evolution predicts, young pulsars can be found orbiting around 
high-mass companion stars (e.g., van den Heuvel & Rappaport, 1987). The discussion 
of Sects. 3 and 4 applies to these systems and time variable high energy emission can 
be produced. The ratio r^jRic is in the range 10“^ — 10®, i.e., substantially smaller 
than the Crab Nebula value. It is therefore possible to test a region of parameter 
space in principle different from that one of the Crab Nebula. 

The prototypical example of a powerful pulsar in a high-mass system is the 
PSR 1259-63 system (Johnston et al.^ 1992). The PSR B1259—63 system provides 
an ideal laboratory to study the interaction of pulsar winds and circumbinary mate¬ 
rial from companion stars. PSR B1259—63 is the hrst radio pulsar found orbiting a 
massive non-degenerate stellar companion. The pulsar has spin period P = 47.76 ms, 
period derivative P = 2.279 • 10“®®, and spindown luminosity Er ~ 9 • 10®® erg s“® . 
The orbital period is Porb = 3.4 yrs, and the orbital eccentricity e = 0.87. Its compan¬ 
ion is SS 2883, a 10th magnitude Be star of luminosity Lo = 5.8 • 10® with radius 
R* ~ 11 i?@ 8 • 10®® cm and an estimated distance from Earth of 2.5 kpc. Be stars 

are characterized by high density mass outflows in the stellar equatorial plane and 
the wide orbit and high eccentricity of the PSR B1259—63 system provide a highly 
variable environment for the interaction between the pulsar and Be star winds. Ac¬ 
cretion near periastron is unlikely for a broad range of reasonable outflow parameters 
(T,n') of the PSR B1259—63 system (TAK94). 

Preliminary results of a multiwavelength campaign organized to monitor the pe¬ 
riastron passage of PSR 1259-63 (January 9, 1994) are consistent with time variable 
pulsar-driven shock emission in the system. The radio pulsar disappeared around De¬ 
cember 21, 1993 and reappeared around February 4, 1994, indicating the presence of 
absorbing/dispersing circumstellar material near the periastron region (Manchester, 
1994). Three ASCA observations of this system in the energy range 0.5-10 keV re¬ 
vealed an unpulsed source with variable intensity of order of 10®® erg s“® and variable 
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spectrum with photon index a ~ 1.7 — 2 (Tavani eA al., 1994a; Kaspi eA al., 1994). 
The only CGRO instrument capable of detecting a faint hard X-ray source (at the 
level of a few mCrab) is OSSE. A very signihcant emission was detected by OSSE 
from PSR B1259—63 in the energy range 30-300 keV during a 3-week observation® 
covering the periastron passage (Grove et al., 1994). A preliminary analysis gives 
an OSSE spectrum consistent with the power-law emission detected by ASCA. The 
characteristics of the X-ray and hard X-ray emission of PSR B1259—63 appear to be 
consistent with those expected from shock emission. The lack of pulsations of the 
X-ray emission (checked independently by ASCA, BATSE and OSSE) argues against 
accretion. Only nebular shock emission or emission due to a ‘propeller effect’ can 
explain the observations (Tavani et al, 1994b). The hardness of the spectrum and 
its power-law nature argue in favor of nebular shock emission. We notice that the 
upper limits established by COMPTEL and EGRET for emission from ~ 1 MeV to 
~ 30 GeV (Tavani et at., 1994b) marginally exclude a shock emission mechanism 
dominated by IC cooling (case B of Sect. 4.1.2). The overall efficiency of conversion 
of pulsar spindown luminosity into broad band emission in the X-ray energy range is 
observed to be 0.01 ~ ~ 0.1. 

PSR B1259—63 is not a very powerful pulsar (it is about 100 time less powerful 
than the Crab pulsar). Furthermore, PSR B1259—63 has a long orbital period, and 
is rarely near its companion star. Systems which are more favorable from the point 
of detecting time variable nebular shock emission may yet be discovered. High-mass 
binaries containing Crab-like pulsars with Porb ~ 10 — 100 days at the distance of 
2-3 kpc would produce a nebular high energy flux between 10 and 100 times larger 
than what observed in the PSR B1259—63 system. 

5.2 RECYCLED PULSARS IN LOW-MASS BINARIES 

A large fraction (~ 50%) of millisecond pulsars are in low-mass binaries with white 
dwarf or lower main sequence companions (e.g., Batthacharya & van den Heuvel, 
1991). This is in agreement with what expected from recycled old pulsars spun up 
during an accreting phase in low-mass X-ray binaries (e.g., Ruderman et at., 1989). 
Even in these binaries, high energy nebular emission can be produced because of 
pulsar wind interaction with mass outflows from irradiated companion stars (AT93). 
The best example of a self-stimulated system (with shock emission providing the 
irradiating flux necessary to drive a mass outflow from the irradiated companion, 
Ruderman et at., 1989) is the eclipsing millisecond pulsar PSR 1957+20 (Fruchter 
et at., 1988). Some of these low-mass systems might contain enshrouded millisecond 
pulsars (Tavani, 1991) and high energy shock emission may be the only way to reveal 
their existence. We notice that in these systems, VslRic ~ 10® — 10®, and from Eq. (4) 
the radius Rq can be of the same order as the orbital distance. 

® OSSE did not detect any flux from PSR B1259—63 during an observation carried out in Septem¬ 
ber 1991, i.e., near apastron (Ray et al., 1993). 
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6. The Nature of Galactic Gamma-Ray Sources 

The unidentified 7 -ray sources in the Galaxy may be isolated compact objects 
(such as 7 -ray pulsars, e.g., Thompson, 1994) and binary systems. The 7 -ray luminos¬ 
ity of isolated pulsars is typically in the range 10^^ — 10^® erg s“^ . We notice that the 
relatively large 7 -ray luminosity of the unidentified 7 -sources (with ~ 100 MeV, 

~ 10^® — 10^® erg s“^ for a distance of 2-3 kpc, e.g., Kanbach, 1994) is similar to 
what is expected from pulsar-driven binaries. These unidentified 7 -ray sources do not 
appear to be associated with strong accreting X-ray sources. Furthermore, several of 
them appear to be time variable by a factor ~ 3 within a timescale of weeks-months 
(Fichtel et al., 1994). We argued here that accreting sources are not likely candidates 
for 7 -ray emission. Rather, binaries containing rapidly rotating pulsars are natural 
candidates for high energy emission. A total of about 10 high-mass and low-mass 
systems with rapidly rotating pulsars are predicted within 2-3 kpc from Farth (e.g., 
van den Heuvel & Rappaport, 1987). The high energy shock emission is likely to be 
time dependent and to follow the geometry variation of the pulsar cavity as a func¬ 
tion of orbital phase. Typical orbital periods range from tens to hundreds of days for 
high mass systems and from a few hours to days for low-mass systems. Monitoring 
galactic 7 -ray sources is therefore of crucial importance for their interpretation. 

At present, there is no unambiguous identification of a 7 -ray source showing a 
modulation associated with an orbital period. The hunt for a binary signature of 
7 -ray sources is an important task for the CGRO instruments. An example of the 
work necessary to reveal the nature of 7 -ray sources is the investigation of the galactic 
source 2CG 135+1. Recent EGRET observations reveal that the new 7 -ray error box 
of this source is consistent with the position of the radio-loud massive star LSI 61 303 
which has a distinctive periodicity of the radio continuum emission of 26.49 days 
(Gregory et al.^ 1979; Taylor & Gregory, 1984). Both the intensity and spectrum of 
2GG 135+1 (e.g., Kanbach, 1994; Hermsen, 1994) are consistent with pulsar-driven 
shock emission in a high-mass system (Tavani et al.^ 1994c; see also Maraschi & 
Treves, 1981 for an earlier model). A series of Gompton observations of this source 
are being carried out during 1994 to reveal a possible modulation of the signal with 
the 26.49'^ period. However, a VLA pulsar search performed at different orbital phases 
gave a preliminary negative result^® (Foster et al.^ 1994). It is clear that we need 
more data to resolve this mystery. 

Another class of interesting 7 -ray sources are dense molecular clouds believed to 
produce 7 -rays because of an enhanced target density to cosmic-ray bombardment 
(e.g., Bertsch et al.^ 1993). If this interpretation is correct, the 7 -ray flux should be 
stable in time. However, a recent analysis of EGRET data shows that the molecular 
cloud p Oph (at an estimated distance d ~ 130 pc and spatially coincident with 

^‘^Note, however, that the VLA pulsar search for LSI 61 303 was not sensitive to a highly-dispersed 
Crab-like pulsar (Foster et al., 1994). 
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the 7 -ray source 2CG 353+16) contains two time variable point-like 7 -ray sources in 
addition to diffuse emission (Hunter et al.^ 1994). A radio-loud quasar (PKS 1622- 
253) might be the responsible for one of the two sources (Hunter et al.^ 1994). 
However, the p Oph cloud is a site of star formation (e.g., Wilking, 1985) and young 
binary systems (either containing massive stars with colliding outflows or rapidly 
rotating pulsars) in the Ophiucus cloud should also be considered as candidates for 
7 -ray emission. A signiflcant contribution to the local cosmic ray flux may originate 
inside the molecular cloud (because of young pulsars) rather than outside the cloud. 
This may be a general phenomenon, important for the interpretation of enhanced 7 - 
ray emission from molecular clouds. It is therefore important to test for the existence 
of time variable 7 -ray sources in other molecular clouds. Young high-mass systems 
with energetic pulsars are likely to be associated with star forming regions and might 
be discovered as time-variable 7 -ray sources. Gamma-ray sources in molecular clouds 
may be tracers of young pulsars in binaries. A deep pulsar search in p Oph and 
similar clouds is also encouraged. 
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